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ABSTRACT 

Gamma Ray Bursts (GRBs) are a powerful probe of the high redshift Universe. We present a 
tool to estimate the detection rate of high —z GRBs by a generic detector with defined energy 
band and sensitivity. We base this on a population model that reproduces the observed prop¬ 
erties of GRBs detected by Swift, Fermi and CGRO in the hard X-ray and 7 -ray bands. We 
provide the expected cumulative distributions of the flux and fluence of simulated GRBs in 
different energy bands. We show that scintillator detectors, operating at relatively high ener¬ 
gies (e.g. tens of keV to the MeV), can detect only the most luminous GRBs at high redshifts 
due to the link between the peak spectral energy and the luminosity (E pea k — T; so ) of GRBs. 
We show that the best strategy for catching the largest number of high-;: bursts is to go softer 
(e.g. in the soft X-ray band) but with a very high sensitivity. For instance, an imaging soft 
X-ray detector operating in the 0.2-5 keV energy band reaching a sensitivity, corresponding 
to a fluence, of ~ 10 -8 erg cm -2 is expected to detect « 40 GRBs yr _1 sr -1 at z ^ 5 (« 3 
GRBs yr “ 1 sr -1 at z fz- 10). Once high-;: GRBs are detected the principal issue is to secure 
their redshift. To this aim we estimate their NIR afterglow flux at relatively early times and 
evaluate the effectiveness of following them up and construct usable samples of events with 
any forthcoming GRB mission dedicated to explore the high -2 Universe. 

Key words: Gamma-ray: bursts — Cosmology: observations, dark ages, re-ionization, first 
stars 


1 INTRODUCTION 

The study of the Universe before and during the epoch of reioniza¬ 
tion represents one of the major themes for the next generation of 
space and ground-based observational facilities. Many questions 
about the first phases of structure formation in the early Universe 
are still open: When and how did first stars/galaxies form? What 
are their properties? When and how fast was the Universe enriched 
with metals? How did reionization proceed? 

Thanks to their brightness. Gamma Ray Bursts (GRBs) can be 
detected up to very high redshift, as already shown by the detection 
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of GRB 090423 at 2 = 8.2 (Salvaterra et al. 2009; Tanvir et al. 
2009) and GRB 090429 at 2 ~ 9.4 (Cucchiara et al. 2011). Thus, 
GRBs represent a complementary, and to some extent unique, tool 
to study the early Universe (see e.g. McQuinn et al. 2009; Am¬ 
ati et al. 2013). A statistical sample of high -2 GRBs can provide 
fundamental information about: the number density and properties 
of low-mass galaxies (Tanvir et al. 2012; Basa et al. 2012; Sal¬ 
vaterra et al. 2013; Berger et al. 2014); the neutral-hydrogen frac¬ 
tion (Gallerani et al. 2008; Nagamine et al. 2008; McQuinn et al 
2008; Robertson & Ellis 2012); the escape fraction of UV pho¬ 
tons from high -2 galaxies (Chen et al. 2007; Fynbo et al. 2009); 
the early cosmic magnetic fields (Takahashi et al. 2011), and the 
level of the local inter-galactic radiation field (Inoue et al. 2010). 


© 0000 RAS 



2 Ghirlanda et al. 


Moreover, GRBs will also allow us to measure independently the 
cosmic star-formation rate, even above the limits of current and fu¬ 
ture galaxy surveys (Ishida et al. 2011), provided that we are able to 
quantify the possible biases present in the observed GRB redshift 
distribution (Salvaterra et al. 2012; Trenti et al. 2013; Vergani et 
al. 2014). Even the detection of a single GRB at z > 10 would be 
extremely useful to constrain dark matter models (de Souza et al. 
2013) and primordial non-Gaussianities (Maio et al. 2012). Finally, 
it has been argued that the first, metal-free stars (the so-called 
PopIII stars) can result in powerful GRBs (Meszaros & Rees 2010; 
Komisarov & Barkov 2010; Suwa & Yoka 2011; Piro et al. 2014). 
Thus GRBs offer a powerful route to directly identify such elusive 
objects and study the galaxies in which they are hosted (Campisi 
et al. 2011). Even indirectly, the role of PopIII stars in enriching 
the first galaxies with metals can be studied by looking to the ab¬ 
sorption features of PopII GRBs blowing out in a medium enriched 
by the first PopIII supernovae (Wang et al. 2012; Ma et al. 2014 in 
prep). 

Specific predictions on the detection rate of high-z GRBs 
(Bromm & Loeb 2002; Gorosabel et al. 2004; Daigne, Rossi & 
Mochkovitch 2006; Salvaterra & Chincarini 2007; Salvaterra et al. 
2008; Butler et al. 2010; Qin et al. 2010; Salvaterra et al. 2012; Lit¬ 
tlejohns et al. 2013) are of fundamental importance for the design 
of future missions. In particular, it is often believed that the pop¬ 
ulation of high redshift GRBs could be easily accessed with soft 
X-ray instruments, due to the (1 + z) _1 cosmological redshift of 
their characteristic peak spectral energy. However, the peak energy 
is strongly correlated with the burst luminosity/energy (Yonetoku 
et al. 2004; Amati et al. 2002) and, therefore, the trigger energy 
band and the detector sensitivity are strongly interlaced. 

The typical prompt emission peak energy 77°]^ of GRBs de¬ 
tected by past/present instruments extends from a few keV up to a 
few MeV. The distribution of the observer frame peak energy ver¬ 
sus redshift of the 185 GRBs with measured spectral parameters 
and redshift z (up to September 2014) is shown in Fig. 1 (see also 
Gruber et al. 2011). The spectroscopically confirmed highest red¬ 
shift GRB 090423 (z = 8.2) has an observer frame -Epeak~90 
keV (von Klienin et al. 2009) as measured by Fermi- GBM. Fig. 1 
shows the absence of GRBs with low measured E°^ k at high red- 
shifts. This is a consequence of the spectral-luminosity correlation 
Speak — iiso (where Li BO is the isotropic equivalent peak luminos¬ 
ity - Yonetoku et al. 2004) and the limited sensitivity of current 
GRB detectors. In other words, we can access, at high redshifts, 
only to the most luminous GRBs that have the highest f7 pea k. This 
compensates the (1 + z)~ 4 decrease of S pea k expected if S pea k 
were independent of luminosity. 

This paper aims to study how these two instrumental prop¬ 
erties can be optimized for detecting the largest number of high 
redshift bursts. We base that on a population synthesis code (§2) 
anchored to the observed properties (prompt emission flux and flu- 
ence and optical afterglow flux) of the GRBs detected by Swift, 
Fermi. BATSE (§3). In §4 we provide a general tool (§4) that allows 
to predict the detection rate of GRBs for a detector with specified 
energy band and sensitivity. Beside their detection, we also study 
the possibility to follow up the high-z GRB population in the X- 
ray and NIR bands through the characterization of their afterglow 
flux level (§4). A standard cosmology with h = Ha = 0.7 and 
Hm = 0.3 is assumed. 
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Figure 1. Observer frame peak energy of the vF v spectrum versus redshift 
of 185 GRBs with measured redshift z and well determined prompt emis¬ 
sion -Epeak- Different symbols (as displayed in the legend) correspond to 
detectors/instruments that provided the measure of E pea k. 


2 SIMULATION SETUP 


The main assumptions to simulate the burst population are: 


(1) GRBs are distributed in redshift (up to z = 20) following 
the GRB formation rate: 

Vt(z) oc (1 + z) 5 t/>*(z) (1) 


where t/>*(z) is the comoving cosmic star formation rate (Li 2008; 
see also Hopkins & Beacom 2008): 


0.0157+ 0.118z 
^ Z) ~ 1 + (z/3.23) 4 ' 66 


( 2 ) 


where t/t(z) has units of Mg yr -1 Mpc -3 . Salvaterra et al. (2012 
- S12 hereafter) have shown that in order to reproduce the observed 
z distribution of the Swift complete sample and the BATSE count 
distribution 5 = 1.7 is needed. We adopted this value. This model 
does not include the possible contribution at z ^ 10 of GRBs pro¬ 
duced from popIII stars. The density evolution of the GRB forma¬ 
tion rate is supported by two considerations: (a) it is consistent with 
the metallicity bias which appears from recent studies of the GRB 
hosts (e.g. Vergani et al. 2014), and (b) it is consistent with the 
expectation of the collapsar model. 

(2) The assumed GRB luminosity function (LF) is a broken 
power law (as derived in SI2): 


Ljso f L cut 

L iso > L cut (3) 

[a, b, Lcut] = [-1.50,-2.32,3.8 x 
10 52 ergs^ 1 ] were constrained by S12, through the BAT6 sample. 
The LF is assumed to extend between minimum and maximum lu¬ 
minosities 10 49 erg s _1 and 10 55 erg s _1 (Pescalli et al. 2014). 
The results for the rate of high-z GRBs are almost insensitive to 
the choice of these bracketing values since the LF is steeply decay¬ 
ing at high luminosities and the low luminosity limit affects only 
the dim bursts that are most likely unseen by any conceivable de¬ 
tectors. 

An alternative model (e.g. Salvaterra et al. 2009; 2012) assumes 
that what evolves with redshift is the break L cu t of the luminosity 


^(I-nso) 


Liso 

Lcut 

Liso 


. Leu t , 

where the parameter values 
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function rather than the GRB formation efficiency. We have tested 
also this assumption and show in §5 that our conclusions are un¬ 
changed. 

(3) To compute the flux and fluence of GRBs in a given energy 
range, we need to associate a spectrum to each simulated burst. As¬ 
suming that all bursts have a spectrum described by a double power 
law function (Band et al. 1993) smoothly joined at the peak energy 
Ep ea k, we assign to each simulated burst a low energy photon spec¬ 
tral index a and a high energy photon spectral index /3 drawn from 
Gaussian distributions centered at —1 and —2.3 with a = 0.2, re¬ 
spectively (e.g. Kaneko et al. 2006; Nava et al. 2011; Goldstein et 
al. 2012). 

(4) The spectral peak energy Speak is obtained through the 
Speak — Si 3 o correlation (Yonetoku et al. 2004). Li ao represents 
the peak isotropic luminosity, i.e. related to the flux at the peak 
of the 7-ray light curve. We adopt the S pea k — Si so correla¬ 
tion obtained with the complete BAT6 sample (Nava et al. 2012): 
log (Speak) = —25.33 + 0.53 log (I/iso). We account for the scat¬ 
ter a = 0.29 dex of the data points around this correlation (Nava et 
al. 2012). 

(5) Similarly, in order to have also a fluence associated to each 
simulated burst, we adopt the S pea k — Si so correlation (Amati et 
al. 2002) between the peak energy and the isotropic equivalent en¬ 
ergy Si 30 , as derived with the BAT6 sample (Nava et al. 2012): 
log (Speak) = —29.6 + 0.61 log (Si so ). Also in this case we con¬ 
sider the scatter a = 0.25 dex. 

Note that our simulation approach does not need to specify the 
distribution of S pea k front where assign the values of this param¬ 
eter to the simulated bursts. Indeed, the starting assumption of 
the simulation is the luminosity function. Therefore, the simulated 
GRB population will have a distribution of Li ao which is the form 
of the assumed luminosity function. It is the assumed Speak — L lso 
correlation which transforms the Li ao distribution into the distribu¬ 
tion of Speak. 

The latter passage is not a 1:1 mapping since we allow for the 
scatter of the correlations. For both the S pea k — Si so and the 
Speak — Li ao correlation the scatter that we adopt represents the a 
of the gaussian distribution of the distance of the data points mea¬ 
sured perpendicular to the best fit line of the correlations. It is still 
debated if and how much these correlations are affected by selec¬ 
tion effects (Band & Preece 2005; Nakar & Piran 2005; Butler et 
al. 2007; Butler, Kocevski & Bloom 2009; Shahmoradi & Nemiroff 
2011; Kocevski 2012) despite some tests suggest that this might 
not be the case (Bosnjak et al. 2008; Ghirlanda et al. 2008; Nava 
et al. 2008; Amati, Frontera & Guidorzi 2009; Krimm et al. 2009; 
Ghirlanda et al. 2012; Dainotti et al. 2014; Mochkovitch & Nava 
2014). However, there seems to be an agreement that, if not prop¬ 
erly a correlation, there could be a boundary in the Speak — Ei so 
and Speak — L^o plane which divides the plane in two regions: 
the right hand side of the correlations lacking of bursts with large 
Siso and Li S o (whose lack might be due to some physical reason - 
e.g. Nakar & Piran 2005) and the left hand side of the correlations 
(corresponding to GRBs with low Si ao and Li ao ) which might be 
affected by instrumental biases. If so, the plane could be uniformly 
filled with bursts on the left hand side of the correlations (i.e. at in¬ 
termediate/low Siso, Li ao ) with no correlation between Si ao and 
Speak and Li a o and S pea k. A fraction of these events could be 
missed due to some instrumental selection effects that shapes its 
cut in the plane so to make an apparent correlation. We have tested 
also this hypothesis and discuss it in §5. 


With these assumptions we derive the 1-s peak flux P of each 
simulated burst in a given energy range A E = [Ei, E2]: 


Liso iff N{E)dE 

47 Td ^z) 2 ' f£™H]+*>EN{E)dE 


where N(E) is the observer frame photon spectrum and cIl(z) is 
the luminosity distance corresponding to the redshift 2. The fluence 
F is: 


Eiso(l + z) Sei E N(E)dE 

4nd L {z)* ' J^ i f+^EN(E)dE 


The integrals at the denominators are performed over the rest 
frame IkeV-lOMeV energy range, which is typically adopted to 
derive Ei ao and Li ao . 


3 SIMULATION NORMALIZATION AND 
CONSISTENCY CHECKS 

The simulated population of GRBs is normalized to the actual pop¬ 
ulation of bright bursts detected by Swift/BAT. The BAT6 sample 
of S12 considered GRBs detected by Swift! BAT with peak flux 
P ^ 2.6 ph cm -2 s -1 (integrated in the 15-150 keV energy range) 
selected only for their favorable observing conditions (Jakobsson et 
al. 2006). The sample contains 58 GRBs and is 95% complete in 
redshift. 

We have extracted from the current Swift sample (773 events 
detected until July 2014 with Tgo ^ 2 s) 1 the 204 GRBs with 1- 
s peak flux P ^ 2.6 ph cm~ 2 s -1 (i.e. the same flux threshold 
used to define the BAT6). Considering a typical field of view of 1.4 
sr and a mission lifetime of approximately 9.5 years, we estimate 
a detection rate by Swift of ~15 events yr~ 2 sr^ 1 with peak flux 
P ^2.6 ph cm~ 2 s -1 . We normalize the simulated GRB population 
to this rate. 

The peak flux cumulative rate distribution of the real Swift 
sample of 204 GRBs with P ^ 2.6 ph cm -2 s -1 is shown in the 
top panel of Fig. 2 (green hatched region representing the statistical 
uncertainty on the rate). Swift has detected many more GRBs with 
peak flux smaller than this limit which is about six times larger 
than the actual Swift flux limit. The consistency of the simulated 
population with the entire Swift sample is discussed in §4.1. The 
corresponding distribution of the simulated GRB population is rep¬ 
resented by the cyan solid line in Fig.2. As shown in the bottom 
left panel of Fig.2, the simulated bursts with -Pris-mo] ^ 2.6 ph 
cm -2 s^ 1 (cyan solid line) have a redshift distribution consistent 
with that of the real BAT6 sample (hatched histogram). 

Note that the comparison of the population of simulated bursts 
with the real swift GRB sample shown in Fig.2 (top panels and 
bottom left panel) are consistency checks. Indeed, their scope is to 
show that the assumptions (§2), on which the code is based, are 
correctly implemented and reproduce the population of bursts (i.e. 
Swift bright events) from which they are derived. 

3.1 The Fermi GBM and BATSE samples 

An interesting test is to verify if the simulated GRB population is 
representative also of larger samples of GRBs detected by other 
instruments besides Swift-BAT. 


1 http://swift.gsfc.nasa.gov/docs/swift/archive/grb_table/ 
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Figure 2. Comparison of the simulated GRB population (solid lines) with real samples of bursts (solid filled and hatched histograms). Top left panel: cumulative 
flux distribution of the bright Swift bursts (hatched green region - one a contingency of the distribution), of the CG/?0-BATSE 4B catalog (solid filled grey 
region) and of the Fermi-GBM 2-year catalog (hatched orange region). Swift fluxes are integrated in the the 15-150 keV range, while BATSE and GBM 
fluxes in the 50-300 keV energy range. The blue, black and cyan solid lines show the flux distribution (in the corresponding energy range) of the simulated 
GRB population corresponding to Fermi , BATSE and Swift. Top right panel: cumulative fluence distributions. Same color coding and symbols of the left 
panel. Bottom left panel: redshift distribution of the BAT6 sample (hatched histogram) and of the simulated sample with the same flux limit P ^2.6 ph cm -2 
s —1 (solid cyan line). Bottom right panel: extinction corrected optical afterglow flux (in the R band) at 1 lh of the BAT6 sample (from Melandri et al. 2014) 
compared to the simulated population with the same peak flux cut (solid cyan line). 


The Gamma Burst Monitor (GBM) on board Fermi has de¬ 
tected many bursts and systematic analysis of their spectral prop¬ 
erties exist (Nava et al. 2011; Goldstein et al. 2012; Gruber et al. 
2014). The first two years GBM sample (Goldstein et al. 2012) con¬ 
tains 478 GRBs with the corresponding spectral parameters, peak 
flux P and fluence F. We considered the flux integrated in the 50- 
300 keV energy range and cut the GBM sample at P ^ 0.7 ph 
cm' 2 s” 1 in order to overcome the possible incompleteness of the 
sample at lower fluxes, obtaining 293 GBM bursts. The GBM is an 
all sky monitor that observes on average ~70% of the sky. There¬ 
fore, the average GBM detection rate of bursts, with P > 0.7 ph 
cm' 2 s^ 1 , is 17 events yr -1 sr _1 . Fig. 2 (top left panel) shows 
that the 50-300 keV flux distribution of real GBM bursts (hatched 
orange region) is consistent with that of the simulated population 
(blue solid line). 

We have considered also the BATSE 4B catalog (Meegan et al. 
1998) containing 1496 bursts with estimated peak flux P and flu¬ 
ence F (both integrated in the 50-300 keV energy range). Also in 
this case we cut the sample to P ^ 0.7 ph cm -2 s -1 (i.e. the same 
level of the GBM), to overcome for incompleteness at low fluxes, 
and obtain 964 events. Considering an average ~75% portion of the 


sky observed by BATSE, the detection rate is ~22 events yr -1 sr _1 
above the considered flux limit. The BATSE peak flux distribution 
(shaded grey region in Fig. 2) is consistent with the simulation (blue 
solid line in the same figure). 

The top right panel of Fig. 2 shows the cumulative fluence dis¬ 
tributions of the Swift, BATSE and GBM samples as shaded regions 
(representing the la statistical uncertainties). The peak flux and flu¬ 
ence distributions of GBM and BATSE are similar (see Nava et al. 
2011a) 2 . The solid lines in the top right panel of Fig. 2 represent 
the simulated GRB population. 


3.2 The afterglow 

The study of high redshift GRBs relies on the capabilities of secur¬ 
ing a redshift for these events. This depends on the brightness of 


2 The three catalogs contain some bursts with only their peak flux or only 
their fluence reported. This is the reason why the cumulative peak flux and 
fluence distributions saturate at slightly different levels for the same instru¬ 
ment. 


© 0000 RAS, MNRAS 000, 000-000 


























High-z GRBs 5 



10" 4 10‘ 2 10 ° 10 2 

Peak flux [ph/cm 2 s] in the corresponding energy band 


Figure 3. Cumulative peak flux distributions representing the rate, i.e. the number of GRBs yr -1 sr 2 , with peak flux ^ P. Different energy ranges where 
the peak flux P is computed are shown with different colors (as shown in the legend). The typical uncertainty on the predicted rates (for any energy range) 
is shown by the vertical bar (see §5). Top panel: entire population of simulated GRBs. Bottom panel : population of high redshift GRBs: z ^ 5(10) shown 
by the solid (dotted) line. For the soft energy band 0.2-5 keV (shown by the red line) a typical Njj = 2 X 10 21 cm~ 2 has been assumed to account for the 
absorption. 


their afterglow in the NIR band. In order to study and predict the af¬ 
terglow flux of the simulated GRB population, we implemented in 
the code the afterglow emission model as developed in van Eerten 
(2010, 2011). This corresponds to the afterglow emission from the 
forward shock of a decelerating fireball in a constant interstellar 
density medium. We assume a constant efficiency r/ = 0.2 of con¬ 
version of kinetic energy to radiation during the prompt phase and 
a typical jet opening angle of 7° (corresponding to the mean value 
of the reconstructed distribution of (9j e t - Ghirlanda et al. 2013). 
For the ISM density we assume a Gaussian distribution centered at 
n = 3 cm -3 and extending from 1 to 30 cm -3 . 

We use the afterglow R flux distribution of the BAT6 sample 
(Melandri et al. 2014 - corrected for extinction, see also Covino et 
al. 2013) to fix the parameters of the external shock micro-physics 
(e.g. Sironi, Spitkovsky & Arons 2013): the fraction of energy that 
the external shock shares with electrons and magnetic field, t e and 
es and the slope of the shock accelerated electrons p. Assuming 
e e = 0.02, tB = 0.008 and p = 2.3 (which also reproduce the 
radio afterglow flux of the BAT6 sample - Ghirlanda et al. 2013b) 
we show that the R flux (computed at 1 lh) distribution of the sim¬ 
ulated bright bursts (i.e. with P ^ 2.6 ph cm -2 s~ 3 - cyan solid 
line in the bottom right panel of Fig. 2) is consistent with that of the 
BAT6 sample (hatched green histogram). A Kolmogorov-Smirnov 


test results in a probability of 0.05 that the two distributions are 
drawn from the same parent population. Differently from the other 
panels, the aim of the bottom right panel of Fig.2 is that of fix¬ 
ing the free afterglow parameters of the simulated GRB population 
in order to match the observed flux distribution of the BAT6 sam¬ 
ple of GRBs (solid filled histogram). We do not pretend to explore 
the afterglow parameter space here (which would require detailed 
modelling of individual bursts - e.g. Zhang et al. 2014) but rather fix 
the parameters (assumed equal for all simulated bursts) in order to 
predict the typical afterglow flux of the simulated burst population. 
Therefore, we consider the agreement of the simulated (solid line) 
and real GRB population afterglow flux distribution shown in the 
bottom panel of Fig.2 satisfactory. With these parameters we can 
predict the afterglow emission of the simulated GRB population at 
any time and any frequency. 


4 RESULTS 

GRBs are detected by current instruments as a significant increase 
of the count rate with respect to the average background. This is 
the rate trigger and it is typically implemented by measuring the 
count rate on several timescales and different energy bands (e.g. 
Lien et al. 2014 for Swift-BAI). Swift-BAI also adopts the image 
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Figure 4. Cumulative fluence distributions representing the rate, i.e. the number of GRBs yr 1 sr —1 , with fluence A F. Different energy ranges where the 
fluence F is computed are shown with different colors (as shown in the legend). Top panel: entire population of simulated GRBs. Bottom panel : population of 
high redshift GRBs: z ^ 5(10) shown by the solid (dotted) line. For the soft energy band 0.2-5 keV (shown by the red line) a typical Ay./ = 2 X 10 21 cm -2 
has been assumed to account for the absorption. 


trigger mode: this searches for any new source in images acquired 
with different integration times. In general both triggers are satis¬ 
fied by most GRBs but the image trigger is particularly sensitive 
to dim/long GRBs: in the Swift GRB sample the relatively small 
fraction of bursts with only image triggers have a smaller peak flux 
than the rest of the bursts (Toma et al. 2011). The time dilation ef¬ 
fect at particularly high redshifts makes GRB pulses appear longer. 
Therefore, image trigger can better detect high-z events. 

The rate trigger is particularly sensitive to the “spikiness” of 
GRBs and, therefore, to the peak flux P, while the image trigger 
is sensitive to the time integrated flux of the burst, i.e. its fluence. 
For these reasons, we present here the peak flux P and fluence F 
distributions of the simulated GRB population. 

Both P and F can be computed on a specific energy range 
A E in our code (Eq. 4, 5). To present our results, based on the 
past/present GRB detectors which were mainly either collimated 
scintillators or imaging coded masks, we have considered the 50- 
300 keV energy range where both BATSE and GBM operate (repre¬ 
sentative of scintillator detectors) and the 15-150 keV energy range 
of Swift-BACT (representative of coded masks). For the predictions 
of the detection of high-z GRBs we considered, for reference, the 
soft energy range 2-50 keV (the BeppoSAX Wide Field Cameras - 
Jager et al. 1997 - were actually operated in the 2-20 keV energy 


range and similarly the Wide X-ray Monitor - Kawai et al. 1999 - 
on board Hete-II), where a coded mask instrument could operate 
and the very soft range 0.2-5 keV, appropriate for an instrument 
like Lobster (Osborne et al. 2013). 


4.1 Peak flux and Fluence 

We have computed the peak flux in the four energy bands given 
above for the entire population of simulated GRBs. The cumula¬ 
tive peak flux distributions are shown in Fig. 3. We show the rate 
of GRB in units of yr -1 sr _1 . The total rate (where the curves sat¬ 
urates at low fluxes) is ~ 1600 GRBs yr -1 sr _ . The flux where 
the curves saturates depends on the energy range where P is com¬ 
puted. The curve showing the flux in the 0.2-5 keV band (red line) 
is slightly different with respect to the other bands. Indeed, to ac¬ 
count for the absorption in the soft X-ray band, we have applied 
a systematic correction of the flux by a factor 1.32, which corre¬ 
sponds to the ratio between the intrinsic and absorbed flux for a 
typical power law photon spectrum ( N(E ) oc E~ v ) with photon 
index T = 1 and assuming an average Nh = 2 x 10 21 cm~ 2 
at z = 0 (Evans et al. 2009; Willingale et al. 2013). The bottom 
panel of Fig. 3 shows the cumulative peak flux distribution for the 
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population of GRBs at z ^ 5 (solid lines) and at z ^ 10 (dotted 
lines). 

The curves in Fig. 3 can be used to predict the detection rate 
of an instrument operating in any of the four energy bands con¬ 
sidered once its flux limit is known. A note on the use of Fig. 3: 
since fluxes in different energy ranges are represented, the thresh¬ 
old (i.e. sensitivity limit of a given instrument) should be applied to 
the corresponding curve. 

For instance, assuming that the flux limit of Swift-B AT is ~0.4 
ph cm ~ 2 s -1 in the 15-150 keV energy range, the top panel of 
Fig. 3 shows that (cyan line) this limit corresponds to ~56 GRBs 
yr -1 sr^ 1 . Considering an average field of view of 1.4 sr and a 
9.5 yr mission lifetime, this rate corresponds to 745 GRBs detected 
by Swift with P ^ 0.4 ph cm -2 s -1 (computed in the 15-150 
keV energy band). This is in agreement with the ~730 real GRBs 
(with T 90 ^ 2 s) present in the Swift on line catalog with peak flux 
above this limit. Similarly one can use the bottom panel of Fig. 3 
to compute the rate of GRBs at high redshift that can trigger Swift. 
We predict that the rate of GRBs at 2 ^ 5 detectable by Swift at 
a level of P ^0.4 ph cm -2 s -1 in the 15-150 keV should be 1.2 
yr -1 sr~ x (corresponding to ~ 16 GRBs over the present mission 
lifetime consistent with the observed number). 

In the Fermi GRB population, instead, assuming a flux limit 
of 0.7 ph cm -2 s -1 in the 50-300 keV energy range and consid¬ 
ering the corresponding (blue) curve in Fig. 3, we predict a rate of 
15 yr ^ 1 sr -1 over all redshifts (top panel of Fig. 2), and a rate of 
0.1 yr -1 sr -1 at 2 ^ 5 (bottom panel of Fig. 2). These rates are 
consistent with the Fermi detection rate above this flux limit. 

Fig. 4 shows the fluence cumulative distribution of the sim¬ 
ulated population. Color codes are as in Fig. 3. We note that the 
curves corresponding to different energy ranges are very similar at 
relatively low fluence values. This is because we show here the en¬ 
ergy fluence, i.e. in units of erg cm -2 (as typically reported in GRB 
databases). In the appendix (Fig. Al) the same cumulative curves 
are shown with the photon fluence (i.e. in units of ph cm -2 ) which 
might also be useful in estimating the detection rate of GRBs for an 
instrument working with image triggers. As explained for the peak 
flux, each curve represented in Fig. 3 should be compared with the 
fluence limit in the corresponding energy range. 
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Figure 5. Observer frame peak energy versus redshift. Real GRBs with 
measured z and well constrained 77 pe ak (185 GRBs) are shown with dif¬ 
ferent symbols (as described by the legend). The grey cloud of points rep¬ 
resents the simulated GRB population. The dark(light) grey contour repre¬ 
sents the 68(95)% containment region of the simulated bursts. The plot is 
truncated at z = 9 for clarity. 
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4.2 High energy detector triggers 

Current GRB detectors, operating at energies ^10 keV, are detect¬ 
ing the most luminous GRBs at high redshifts. Indeed, due to the 
Speak — Lise correlation, these events have relatively high peak en¬ 
ergies. Fig. 5 shows the observer frame peak energy of real GRBs 
detected by different instruments (as labeled) versus their redshift 
(same as Fig. 2) together with the entire simulated GRB popula¬ 
tion (grey region). At high 2 the detection of the most luminous 
bursts, due to the S pea k — Li so correlation, also selects the GRBs 
with the largest peak energies. This effect can also be seen in Fig. 6 
which shows the observer frame peak energy versus the bolometric 
flux for the real 185 GRBs with well constrained spectral parame¬ 
ters and for the simulated GRB population (same symbols of Fig. 
5). Here we also highlight (yellow region) the simulated popula¬ 
tion of GRBs at 2 ^ 5. We note that, due to the E pea k — Liao 
correlation, a correlation is present in the observer frame between 
the observed peak energy and the peak flux (e.g. Nava et al. 2008). 
Through Fig. 6 one can understand that actual GRB detectors, most 
sensitive from 10 keV up to a few MeV, select those bursts with 
peak energy in this energy range. These events, due to the exis¬ 
tence of the D P eak — Liao correlation, are also the brightest (and 


Figure 6. Observer frame peak energy versus bolometric flux. Real GRBs 
with measured z and well constrained prompt emission spectral peak energy 
(185 GRBs) are shown with different symbols. The grey cloud of points 
represents the simulated GRB population. Yellow points are GRBs at 2 > 
5. 


most luminous) bursts although their number is shaped by an over¬ 
all decreasing luminosity function. 

Fig. 6 shows that in order to efficiently detect high redshift 
GRBs (yellow symbols) one needs a detector working in a softer 
energy band (e.g. 0.1-10 keV) but also much more sensitive than 
current detectors. 

We have considered four energy bands A E based on two 
classes of detectors: imaging coded mask detectors, typically oper¬ 
ated up to few hundreds of keV and scintillator detectors operating 
in 10 keV-1 MeV range. It is clear that the hunt for high -2 GRBs is 
optimized with detectors operated in the soft energy range, i.e. few 
keV. However, it is worth exploring the rate of detections of GRBs 
(and high -2 bursts) by both types of detectors. Indeed, for a config¬ 
uration with both these detectors, for instance, while the detection 
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Figure 7. Cumulative peak flux (in the 50-300 keV energy range) distri¬ 
bution for different cuts of the soft (0.2-5 keV) energy fluence (as shown 
for the four solid lines). The solid, dashed and dotted lines show the entire, 
z ^ 5 and z ^ 10 population of simulated bursts. The four lines for each 
set correspond to the 0.2-5 keV fluence cuts of F ^ 10 -7 , 10 -8 , 10 -9 
and 10“ 12 erg cm -2 . These curves serve to estimate the rate of common 
triggers for two instruments: one working in the 0.2-5 keV energy range 
and one in the 50-300 keV energy range, once both have assigned fluence 
and peak flux thresholds, respectively. 

rate of high-z bursts is maximized by the soft energy detector, the 
study of the GRB prompt emission properties (spectral shape, lu¬ 
minosity etc.) is possible for those bursts also detected by the high 
energy detector. 

Fig. 7 shows the 50-300 keV peak flux cumulative rate distri¬ 
bution of the simulated GRB population for different cuts of their 
0.2-5 keV energy fluence (as labeled). We also show the rate for 
GRBs at z ^ 5 and z ^ 10 (blue dashed and red dotted lines, 
respectively). For a soft X-ray detector (working in the 0.2-5 keV 
energy range) with a sensitivity threshold of 10 -8 erg cm ' 2 cou¬ 
pled with a typical scintillator detector (working in the 50-300 keV 
energy range) with a sensitivity of 10 -1 phot cm ' 2 s - , the rate 
of GRBs detected by both detectors is ~65 yr -1 sr -1 and 2 yr -1 
sr -1 at z ^ 5. 

4.3 Afterglow 

In addition to the detection of the high redshift GRB population, 
it is fundamental to assess their redshifts. This is possible through 
early time NIR follow up observations. To this aim we have com¬ 
puted the NIR flux of the simulated GRB population as described 
in §3. Fig. 8 shows the cumulative J band flux distribution for the 
GRBs at z ^ 5 (solid lines) and z ^ 10 (dashed lines) considering 
different thresholds of their 0.2-5 keV energy fluence. 

A way to optimize the follow up of high redshift GRBs could 
be with a dedicated NIR telescope on board any GRB dedicated 
satellite trough rapid slewing. Assuming aim telescope in space 
we expect H~22.5 in 500 seconds of integration and H~22 in lh 
for imaging and low (R=20) resolution spectroscopy. Higher reso¬ 
lution spectroscopy (e.g. R~1000), needed for metallicity studies, 
can be performed for afterglows brighter than H~18. With these 
estimates we evaluate, for imaging and assuming a rapid follow up 
i.e. within 500s, a rate of ~40 GRBs yr -1 sr -1 at z ^ 5 and 3 yr -1 
sr -1 at z > 10 (with 0.2-5 keV fluence > 10 -8 erg cm -2 ). Metal¬ 
licity studies could be performed for ~30(1.5) GRBs yr -1 sr -1 at 



Figure 9. Cumulative distributions of the (0.3-10 keV) X-ray flux at 5 h 
after the burst. Same symbols as in Fig. 5. 

z ^ 5(10) (considering a GRB detector with a limiting 0.2-5 keV 
fluence ^ 10 -8 erg cm -2 ). 

Finally we derive the X-ray (0.3-10 keV) flux of the simu¬ 
lated GRB population. This is interesting for the use of GRBs as 
probes of the chemical composition and properties of the interven¬ 
ing medium (e.g. Campana et al. 2011) through the metal absorp¬ 
tion edges imprinted on their X-ray spectra. 

The X-ray properties of the GRBs of the BAT 6 sample have 
been presented in D'Avanzo et al. (2012). They found that the X- 
ray luminosity is correlated with both the prompt emission 7 -ray 
isotropic luminosity and energy (Lj so and E lso respectively - Am- 
ati et al. 2002; Yonetoku et al. 2004). There is evidence that the 
early X-ray emission observed in GRBs has a different origin to 
the afterglow produced at external shocks (Ghisellini et al. 2009). 
Besides there could be a connection between the prompt and X-ray 
emission (Margutti et al., 2013). In particular, the early time emis¬ 
sion in the X-ray band could be due to an emission component 
related to the long lived activity of the central engine (invoked also 
by late time X-ray flares observed in the Swift XRT light curves 
- Burrows et al. 2007). Therefore, we estimated the X-ray flux of 
the simulated GRB population through the empirical correlations 
found by D’Avanzo et al. (2012). Fig. 9 shows the cumulative dis¬ 
tribution of the flux in the 0.3-10 keV energy range for different 
cuts of the 0.2-5 keV fluence. 

With a Target Of Opportunity (TOO) reaction time of 4-5 
hours Athena (Nandra et al. 2013) will be able to provide high res¬ 
olution spectroscopy for a large number of GRB afterglows. In or¬ 
der to enable high resolution spectroscopy in the microcalorimeter 
(XIFU) instrument, that delivers 2.5 eV resolution, a total fluence 
of 10 5 cts is needed (lonker et al. 2013). This corresponds to typi¬ 
cal afterglow flux of ^ 10 -12 erg cm -2 s -1 with a integration time 
of ~50-100 ksec. With this flux limit, the expected rate of GRBs is 
3(0.2) yr -1 sr -1 at z ^ 5(10). 

The effective number of GRB that can be followed up by 
Athena wil depend on the field of view and the observational strat¬ 
egy of the GRB mission dedicated to detect GRBs and provide 
trigger alerts for the X-ray observations. For a limited (^1-3 sr) 
field of view of a soft X-ray detector (operated e.g. in the 0.2-5 
keV), the optimal strategy to maximize the number of follow-up 
with Athena will be to point towards the sky accessible by Athena 
at that given time (taking into account its observing program). This 
strategy would mimimize the slewing time in case of a trigger. In 
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Figure 8. Cumulative NIR flux distribution of the afterglow of GRBs at z > 5 and z > 10 (blue and red lines respectively) at 500 s for different cuts on the 
0.2-5 keV fluence. For 5 < z < 10 we considered the J band and for z > 10 the H band. We assume Av = 0 for the population of high redshift GRBs. The 
dotted lines show the J flux at 1 h. 


such a case the Athena observing efficiency will be close to 100% 
and the TOO reaction time would be also minimized. For a much 
larger FOV, (> 1/2 sky), the number of GRB follow up by Athena 
will, on the contrary reduce to 50%. 


5 DISCUSSION 

We have presented the prompt and afterglow properties of a syn¬ 
thetic population of long GRBs. The main assumptions of the sim¬ 
ulations are the luminosity function and the GRB formation rate 
(proportional to the star formation rate). GRB prompt emission 
properties are derived through the empirical correlations between 
the peak energy of their spectra and the isotropic equivalent energy 
and luminosity (E pea k — Ei ao and E pea k — Lj so correlations). The 
simulated population is normalized to the bright sample of GRBs 
detected by Swift (i.e GRBs with peak flux Eis-isokeV ^ 2.6 ph 
cm' 2 s _1 ). The simulated population also reproduces the proper¬ 
ties (redshift distribution and R afterglow flux distribution) of these 
real bursts. We have shown that the simulated population is repre¬ 
sentative also of the larger population of bursts detected by Fermi 
and BATSE. 


The detection rates presented are based on the best fit values 
of the assumed GRB formation rate (Eq. 2) and their luminosity 
function (Eq. 3). The main source of uncertainty for the predictions 
of the rates of high redshift bursts is the strength 8 of the evolution 
with redshift of the GRB formation rate. Considering 5 = 1.7±0.5 
(S12), the rates can be larger/smaller by a factor ~1.5. This only 
affects the absolute rate of detection and not the comparison of the 
fraction of bursts accessible with instruments working in different 
energy bands as those presented here (§4.2. 4.3). 

In general our results show that the largest rate of high-z 
GRBs can be detected with a soft X-ray detector (e.g. operated 
in the 0.2-5 keV energy range like a Lobster-type wide-held fo¬ 
cussing X-ray) because this is where the bulk of the high-z GRBs 
have their peak energy (due to the cosmological redshift). FIow- 
ever, such an instrument should be designed so to have a sensitiv¬ 
ity that accounts properly for the the presence of a link between 
the .Epeak and the luminosity/energy of GRBs (E pea k — Li ao and 
Epeak — Eiao correlations). To asses a population of GRBs say a 
factor of ten softer in E pea k requires a sensitivity improvement of 
a factor ~ 100. This is due to the slope of the observer frame distri¬ 
bution of the points in Fig.6 which is a result of the E pea k — Ei ao 
correlation. The purpose of this work is to provide a tool for de- 
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signing instruments for detecting a large sample of GRBs to study 
the high —z Universe. 

The results presented so far are based on the set of assump¬ 
tions described in §2: the evolution of the GRB formation rate with 
redshift (Eq. 1) and the E peak - E iso and E pea , k - L iso corre¬ 
lations which are used, considering also their scatter, to link the 
energetic/luminosity of the simulated bursts with their peak energy. 

The first assumption is motivated by the fact that an evolu¬ 
tion of the GRB density is also consistent with the metallicity bias 
as found from the study of GRB hosts (e.g. Vergani et al. 2014), 
and it is theoretically expected within the collapsar model. How¬ 
ever, also the alternative scenario in which the luminosity function 
evolves with redshift has been shown to be consistent with the ob¬ 
servations (e.g. Salvaterra et al. 2012). For completeness we have 
also explored this possibility, assuming the GRB rate without any 
redshift evolution (i.e. <5 = 0 in Eq.l) but, instead, a luminosity 
function described by a broken powerlaw with a break luminos¬ 
ity L cu t oc (1 + z) 2 ' 1 . This model increases the fraction of lumi¬ 
nous bursts at high redshifts: the results are shown in Fig. 12 by the 
dashed lines. Compared with the results obtained with an evolving 
GRB rate (§2 - shown by the solid lines in Fig. 12), the global ef¬ 
fect is to reduce the rate of dim bursts due to the increase of the 
high redshift and bright bursts (bottom panel of Fig. 12). However, 
these differences are consistent with the typical uncertainty on the 
estimated rates (shown by the vertical bar in Fig. 12). 

The other assumptions, i.e. the Speak — Si so and E peak — Li ao 
correlations can also be tested. We have assumed (§2) these two 
correlations with their scatter. However, it is still debated if they are 
due to instrumental selection effects which, if present, should limit 
the detection of GRBs of intermediate/low luminosity/energy, i.e. 
those on the left hand side of these correlations (in the S pea k — Ei so 
and Speak — Sso correlations we typically plot the x-axis rep¬ 
resenting the isotropic energy Si ao or luminosity Li ao ). Although 
excluded by recent simulations (Ghirlanda et al. 2012) we have 
considered the extreme case of a uniform distribution of simulated 
bursts on the left hand side of the S pea k — Si ao and S pea k — Liso 
correlations. In practice we assume that the S pea k — S so and 
S P eak — Si a o correlations are boundaries in their respective planes 
(e.g. Nakar & Piran 2005; Butler Bloom & Poznanski 2010) which 
are. therefore divided in two regions: on the right hand side of the 
correlations there are no bursts (if they were present we should have 
detected) while on the left hand side there is no intrinsic correla¬ 
tion between e.g. S pea k and Z/i ao (Si ao ). The results obtained with 
this assumption are shown by the dashed lines in Fig. 13. From the 
comparison with the results (solid lines) obtained with the set of 
assumptions of §2 (i.e. assuming instead the correlations with their 
scatter), it appears that the difference is minimal and consistent 
within the uncertainty on the estimated rates. 


6 CONCLUSION 

Through our code we can derive the prompt emission flux and flu- 
ence and the afterglow flux (at any time and frequency) of the pop¬ 
ulation of GRBs distributed up to 2 = 20. We note that our code 
does not include the possibility that the GRB rate at high redshift 
is enhanced by the primordial massive (popIII) stars which could 
even be dominating at 2 ^ 10 (Campisi et al. 2011). We have de¬ 
rived for the first time also the afterglow J and H band flux for high 
redshift GRBs. This is important because it can be used to assess 
the fraction of detected GRBs whose afterglow can be observed (or 


spectrum acquired) with an IR telescope (in order to maximize the 
study of the high -2 population). 

The scope of this paper is to provide a tool to estimate the de¬ 
tection rate of GRBs and, in particular, of high -2 bursts. High -2 
GRBs are of particular relevance for planning forthcoming dedi¬ 
cated missions to study the high -2 Universe. We have shown here 
the flux and fluence rate distributions considering four possible en¬ 
ergy bands where present and future GRB detectors can operate. 
The presented curves can be used to estimate approximately the 
detection rates expected for an instrument operating in a given en¬ 
ergy range and with a given sensitivity threshold. 
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Figure 10. Cumulative distributions of the peak flux (in units of erg cm 2 
s —1 ). Same lines and color coding of Fig.4. 

APPENDIX 

We show in Fig. 10 the cumulative distributions of the peak flux in 
units of erg cm -2 s _1 and in Fig.11 of the fluence in units of ph 
cm -2 . In the results presented in the paper we have assumed a den¬ 
sity evolution of the GRB formation rate with a constant luminosity 
function, i.e. a broken powerlaw (§2). If the evolution is instead as¬ 
signed to the luminosity function, we adopt the model parameters 
of Salvaterra et al. (2012) where there is an evolution of the cut of 
the luminosity function L cu t oc (1 + z) 2 ' 1 . This evolution increases 
the number of high luminosity GRBs at high redshifts. Fig. 12 com¬ 
pares the cumulative peak flux distribution obtained with the as¬ 
sumption of a density evolution of GRBs (§2 - solid lines) and with 
the assumption of a luminosity evolution (dashed lines in Fig. 12). 

The other main hypothesis of our paper is the existence of 
the Epeak — Ei BO and _E pea k — L{ so correlations. It has been dis¬ 
cussed that these could boundaries in the corresponding planes and 
some instrumental/observational effects prevent the exploration of 
the part of the planes which would be uniformly filled with in¬ 
termediate/low Ei so (Li so ) bursts. We have tested this possibility 
with our code and show in Fig.13 our results. We have assumed 
the Speak — Ej so and S pea k — Disc, correlations, as obtained from 
the current data, and assumed that they divide the plane into two 
regions: the right hand side of these correlation is free of bursts, 
since if any GRB would be in this semi-plane it should have been 
detected, while the left hand side of the correlation is filled (uni¬ 
formly) with bursts which have their peak energy uncorrelated with 
the Ei so or Li so . The boundary of the uniformly filled semi-plane 
are defined assuming that bursts have S pea kG (0.01; 10 4 ) keV. 
This different assumption produce a flux cumulative distribution 
(dashed line in Fig. 13) which is only slightly different from that 
obtained assuming the correlations as they are observed with their 
actual scatter (solid lines in Fig. 13). 



Figure 11. Cumulative distributions of the fluence (in units of ph cm 2 ). 
Same lines and color coding of Fig.8. 



Figure 12. Cumulative distributions of the peak flux comparing the results 
obtained assuming the evolution of the GRB formation rate (solid lines - 
§2) and the evolution of the luminosity function (dashed lines). Same lines 
and color coding of Fig.4. 



Figure 13. Cumulative distributions of the peak flux comparing the results 
obtained with the assumptions of §2 (solid lines) and that the -E pea k — -Eiso 
and -Epe a k — T; so correlations are boundaries in the corresponding planes 
(dashed lines). Same lines and color coding of Fig.4. 
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